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ABSTRACT 

We present and discuss the X-ray luminosity function (XLF) of the ROSAT 
Brightest Gluster sample (BGS), an X-ray flux limited sample of clusters of 
galaxies in the northern hemisphere compiled from ROSAT All-Sky Survey data. 

The BGS allows the local cluster XLF to be determined with unprecedented 
accuracy over almost three decades in X-ray luminosity and provides an 
important reference for searches for cluster evolution at higher redshifts. 

We find the significance of evolution in both the XLF amplitude and in the 
characteristic cluster luminosity L*x to be less than l.Scr within the redshift 
range covered by our sample thereby disproving previous claims of strong 
evolution within z ^0.2. 


Subject headings: galaxies: clusters: general — cosmology: observations — 
X-rays: general 


1. Introduction 

The ROSAT Brightest Gluster sample (BGS, Ebeling et ah 1996b, hereafter Paper I) 
is a 90% complete, flux limited sample of the 199 X-ray brightest clusters of galaxies in the 
northern hemisphere {S > 0°), at high Galactic latitudes (|6| > 20°, with redshifts z < 0.3, 
fluxes higher than 4.45 x 10“^^ erg cm“^ s“^ and luminosities higher than 5 x 10^^ erg s“^ in 
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the 0.1-2.4 keV band. Second in size only to the XBACs sample of Ebeling et ah (1996a), 
the BCS is one of the largest statistical clnster samples compiled at X-ray wavelengths to 
date. It is the only large-scale sample available today that is not only X-ray flux limited bnt 
also X-ray selected in the sense that the BCS, nnlike the XBACs, is not limited to systems 
initially fonnd in optical surveys but contains clusters selected by their X-ray properties 
only. 

The BCS thus represents an ideal sample for studies of the formation, distribution 
and evolution of structure on the largest metric and mass scales. Providing important 
constraints on the cosmological parameters governing cluster evolution, the X-ray luminosity 
function (XLF) of clusters of galaxies represents a particularly vital statistic in this context. 

Several previous studies based on much smaller samples of typically 50 clusters found 
the evolution in the cluster X-ray luminosity function to be ‘negative’ in the sense that 
X-ray luminous clusters are more numerous now than they were in the past (Edge et ah 
1990; Gioia et al. 1990; Henry et ah 1992; David et ah 1993). They were, however, not only 
in conflict with other studies, which found no evidence for cluster evolution (e.g., Kowalski 
et al. 1984), but also somewhat inconsistent among themselves. The strong evolution seen 
by Edge et al. (1990) in their sample of 46 X-ray bright clusters at high Galactic latitude 
and < 0.18 is not present in the first two redshift bins (44 clusters at 0.14 < < 0.3) of 

the sample of Gioia et al. (1990) and Henry et al. (1992) who hud significant evolution 
only at z > 0.3. More recently, two studies found no sign of evolution at all in the XLF 
of samples of Abell and AGO clusters at z < 0.36 (Briel & Henry 1993) and < 0.15 
(Burg et al. 1994), respectively. However, these samples were neither X-ray selected nor 
X-ray flux limited and may thus not be fair representations of the cluster population in 
general. At considerably higher redshifts {z > 0.5) on the other hand, studies based on yet 
smaller samples observed in deep X-ray pointings (Bower et al. 1994; Castander et al. 1994; 
Castander et al. 1995) suggest a signihcant drop in the cluster space density as compared 
to the value observed locally. 

Although the overall evidence is thus in favor of negative evolution of the cluster XLF 
at least for X-ray luminous clusters at redshifts well in excess of 0.3, the overall picture is 
anything but clear. With the completion of the BCS we are now able to provide a dehnitive 
answer to the question of whether cluster evolution is signihcant at low to intermediate 
redshifts and, in any case, provide an accurate determination of the local cluster XLF as a 
much-needed reference for ongoing and future evolutionary studies at higher redshifts. The 
implications of our Endings for cosmological models of cluster evolution will be addressed 
in a forthcoming paper (Ebeling et ah, in preparation). We assume Hq = 50 km s“^ Mpc“^ 
and go = 0.5 throughout this paper. 
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2. The BCS XLF and its parametrization 


The BCS as published in Paper I is only 90% complete and corrections for 
incompleteness need to be applied to account for clusters missing from the sample. In 
doing this we use a selection function based on the z, Lx distribution of the serendipitous 
detections in our sample (see §7.1 of Paper I). 


We use the usual dehnition of the unbinned, differential XLF for a sample with flux 
limit fx,hm- 

dn{Lx,z, 


XLF{Lx, z, fx,hm) — 


dL 


X 


where dn{Lx, z, fx,iim) is the space density of clusters with X-ray fluxes above the flux 
limit and X-ray luminosities within an interval dLx around Lx- Since we use an unbinned 
representation, dn is given for each cluster by l/17(Lx, 2 :, /_x,iim), i-e., the inverse search 
volume dehned by the luminosity distance at which the X-ray flux from a cluster with 
intrinsic luminosity Lx would equal the flux limit of the sample. The systems’ X-ray 
temperatures as listed in Paper I are used in the computation of K corrections. 


We use a Schechter function (Schechter 1976) of the form 

^{Lx)=AeM-Lx/Ll)L-^- 

to model the XLF and determine A, L*x and a in a maximum-likelihood £t to the unbinned 
data. 


Table gives an overview of the £t results obtained in all hve standard energy bands 
currently in use within the scientific community: 0.1 — 2.4, 0.5 — 2.0, 0.3 — 3.5, and 2 — 10 
keV, as well as the pseudo-bolometric band from 0.01 to 100 keV. We discuss the results in 
detail in the following paragraphs. 

Figure shows the differential XLF for the BCS in the generic 0.1 — 2.4 keV band 
of the ROSAT observatory. We test the robustness of the £t by comparing it to the XLF 
obtained for the larger, 80% complete BCS (unpublished, cf. Paper I) and find excellent 
agreement. Also shown in Figure are the XLF data points for the high-galactic-latitude 
sample of Edge et ah (1990, hereafter B50) (open diamonds) which is X-ray flux limited 
in the 2 — 10 keV band. Since all 46 B50 clusters have measured X-ray temperatures, the 
conversion of their luminosities to the 0.1 — 2.4 keV band of the BCS is less inaccurate than 
the opposite operation, i.e., the conversion of the BCS luminosities into the 2 — 10 keV 
band, which relies heavily on temperature estimates rather than measured values. To assess 
the impact of band conversion effects we make the comparison between the XLFs of the 
B50 and the BCS in either band. In the 0.1 — 2.4 keV band we End the XLF for the B50 to 
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be, in general, in good agreement with the BCS XLF. At low X-ray luminosities, however, 
the best Schechter function £t for the B50 (the dotted line in Figure and obtained 
with the same ML algorithm used throughout for the fitting of the BCS XLFs) severely 
underpredicts the observed volume density of clusters, indicating incompleteness of the B50 
at Lx ^ 1 X 10"^^ erg s“^ (0.1 — 2.4 keV). This is, however, not surprising given that the 
B50 is, by design, only flux-limited in the 2 — 10 keV band. 

Figure ^ shows the BCS XLF in the ROSAT hard band covering the energy range 
from 0.5 to 2.0 keV. Note that, like for all other energy bands discussed in the following, 
the conversion from the original BCS band (0.1 — 2.4 keV) to the ROSAT hard band was 
performed for each cluster individually using the X-ray temperatures and Galactic column 
densities given in Table 2 of Paper 1 and assuming a metallicity of 0.3 of the solar value. 
Also shown in Figure ^ are the data of the XLF of groups and poor clusters of galaxies 
as presented by Burns et ah (1996, hereafter BLL) and their best power-law XLF. Note 
the good agreement between the two samples as well as the respective XLF models in the 
overlap region between 3 x 10^^ erg s“^ (the lowest luminosity of any BCS cluster) and 
2.6 X 10^^ erg s“^ (the highest luminosity of BLL’s poor cluster sample)]^. Finally, we overlay 
in Fig. H the best fitting Schechter function XLF determined by De Grandi (1996) for the 
111 clusters of the BSGC-KP sample, a subset of a larger ROSAT cluster sample under 
compilation in the southern hemisphere (Guzzo et ah 1995). For luminosities in excess of 
~ 2 X 10^^ erg s“^ De Grandi’s XLF is in very good agreement with the best-fitting BGS 
Schechter function. Below this value the BSGG-KP XLF falls increasingly below both the 
BLL and the BGS fits, reaching a deviation of a factor of 2 and 3, respectively, at the lower 
end of the BSGG-KP luminosity range. 

The reference in the 0.3 — 3.5 keV band is the EMSS cluster sample of Henry et ah 
(1992). Figure ^ shows the BGS XLF in the 0.3 — 3.5 keV band with the EMSS data points 
from the first two redshift shells (0.14 < .^ < 0.2 and 0.2 < z < 0.3) overlaid as shaded 
diamonds. (We remind the reader that the power law descriptions of the EMSS XLF found 
by Henry and coworkers in these shells are consistent with one another within the errors, 
so that the comparison with the whole of the BGS made in Figure § is legitimate.) Note 
the very good agreement between the EMSS and BGS luminosity functions, but also the 
much higher accuracy provided by the BGS. Also shown in Figure ^ is the XLF from the 
third EMSS redshift shell (0.3 < 2 ; < 0.6). Gontrary to our findings for the EMSS XLF for 


^We have tested the compatibility of the two samples by comparing the luminosities of the four clusters 
contained in both the BCS and the BLL sample. We find agreement to within 5 % between the respective 
luminosities for all but one cluster (MKW8) for which the BLL luminosity falls short of the BCS value by 
40%. 
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clusters at 0.14 < z < 0.3, the EMSS XLF of these high-redshift systems differs noticeably 
from the BCS XLF as the local reference. A detailed re-assessment of the signihcance of 
the evolution implied by this discrepancy is presented by Ebeling et al. (in preparation). 

Although, in Table |I|, we do list the results of a Schechter function £t to the BCS data 
in the 2 — 10 keV band, the BCS XLF in this hard energy band should be regarded with 
caution. The 2 — 10 keV band has hardly any overlap with the 0.1 — 2.4 keV BCS detection 
band, which renders a flux conversion that is largely based on estimated temperatures 
a dangerous enterprise. Keeping this caveat in mind we nonetheless End the power law 
description from Piccinotti et al. (1982) to be in good agreement with the BCS XLF 
irrespective of whether the Virgo cluster is included or not (Figure p. 

The best £t parameters from the Schechter function fit to the B50 data published 
by Edge et al. (1990), on the other hand, provide an unacceptable fit to the BCS XLF 
at luminosities in excess of about 1 x 10^® erg s“^, where the £t given by Edge et al. 
falls signihcantly below the BCS data. This failure of the original B50 Schechter function 
£t to describe the BCS XLF is due to an error in the volume calculation in the work of 
Edge and coworkers. Fitting the B50 data with our maximum likelihood algorithm we 
End A = l.SOlgJ® X 10“'^ Mpc“^ (10"^"^ erg s“^)"“^, x 10^“^ erg s~^, and 

a = 1.25lo;2o, in good agreement with the results for the BCS (see Figure § and Table |I|). 


3. Evidence of evolution in the BCS XLF 

We search for evidence of evolution by splitting the BCS at a redshift Zgep and 
independently Etting Schechter functions to the data in the two redshift shells thus created. 
Care has to be taken not to naively misinterpret every statistically signiEcant difierence 
between the best Et parameters in the two shells as signature of evolution. In order to avoid 
efiects due to large scale structure, we vary Zsep from 0.05 to 0.2 and look for trends that are 
robust over a range of Zsep values. Since, due to the Eux-limited nature of our sample, the 
low-luminosity end of the XLF, and thus a, is ill-constrained for the high-redshift subsample 
once ^^sep exceeds 2 ; ~ 0.1, we fix the power law slope a at its overall best-fit value of 1.85 
in the maximum-likelihood fits to the data of either subsample. With a frozen we are thus 
left with two free parameters, the normalization A and the characteristic luminosity L\. 

Figure | shows the contours of the C statistic (which is distributed) of A and L*x 
for some of these low-redshift and high-redshift subsamples of the BCS as a function of Zsep- 
While differences found at Zgep ^ 0.1 can be attributed entirely to large-scale structure, 
a significant decrease in A or at higher values of Zgep would be indicative of negative 
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evolution. As, for 0.05 < .^sep < 0.2, the 68% confidence contours of the low-redshift and 
high-redshift subsamples overlap, we conclude that there is no signihcant evolution in either 
the amplitude A of the cluster XLF or the characteristic luminosity L\ values of Zsep up 
to 0.2. Since, in the high-redshift subsamples with 2;gep ^ 0.16, A becomes ill-constrained 
and increasingly strongly correlated with we also tested for evolution only in L*x by 
holding A constant at its overall best-£t value of 4.33 x 10“^ Mpc“^ (lO^"^ erg s“^)““^, 
a value well within the 68% confidence contours of all hts shown in Fig. We find the 
variations in L*x to be smaller than 30(37)% for 0.1 < Zsep < 0.20(0.22) which is less than 
1.6(1.8 )ct signihcant, conhrming the no-evolution result of Fig. ^ 

As an independent check, we also looked for variations in 17/Knax us u function of both 
z and Lx — and found none. A Kolmogorov-Smirnov test hnds the distribution of K/14iax 
values (whose median is 0.47) to be consistent with uniformity at the greater than 74% 
confidence level suggesting again that the cluster space density of the BCS is homogeneous 
out to the limiting redshift oi z = 0.3. 

4. Conclusions 

Using the ROSAT Brightest Cluster Sample (BCS) as presented by Ebeling et ah 
(1996b) we have established the local X-ray luminosity function (XLF) of clusters of 
galaxies within 2 ; = 0.3 with unprecedented accuracy. We find the XLF to be well 
described by a Schechter function whose free parameters A, L*x, and a we determine in a 
maximum-likelihood fit for all X-ray energy bands currently used within the community. 
Comparing our results with previous measurements of the cluster XLF we find very good 
agreement with the work of Piccinotti et al. (1982), Henry et ah (1992), and Burns et al. 
(1996), as well as with the XLF for the B50 sample of Edge et al. (1990) when the same 
maximum likelihood algorithm is used to determine the best Schechter function £t. 

We hnd no signihcant variations in the amplitude or the characteristic luminosity 
of the best-htting Schechter function as a function of redshift. Also, the distribution of 
U/Umax values is consistent at the 74% conhdence level with a non-evolving space density 
of clusters out to z = 0.3. Our findings do thus not conhrm the claim of strong evolution 
at z ^ 0.2 made by Edge and coworkers but support the notion of Ebeling et al. (1995) 
that the apparent signature of evolution in the B50 sample is due to a combination of its 
high X-ray hux limit in the 2 — 10 keV band and a pronounced, if statistically insignihcant, 
dearth of very X-ray luminous clusters around a redshift of about 0.15. 


H.E. thanks Pat Henry for helpful discussions about maximum-likelihood htting and 
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Table 1. Best fit values of the Schechter function parameters A, and a 


Energy Band 

At 


a 

0.1-2.4 keV 

O.UD_o.46 

9.10+i 

1 oe-+0.09 
-‘-•OO-0.09 

0.5-2.0 keV 

Q 00+0.36 
'^•'^^-0.33 

5.70lJ;i 

-‘-•O'-'-0.09 

0.3-3.5 keV 

4.95b°i« 

10.7b?j 

-1 oo+O.OS 
-*-•^^-0.09 

2-10 keV 

9 qc;+0-22 

12.6b?:2 

1 c;4+0-05 

-‘-•3^-0.06 

bolometric 

f, 41 +0-70 

37.2b^|" 

1 q/I+0.09 

1-o4_o,04 


Tn units of 10 ^ Mpc ^ ( 10 "^^ erg s ^)“ ^ 
bn units of 10^"^ erg s“^ 

Errors given correspond to 68% confidence for 
one interesting parameter (AC* = 1). 
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Fig. 1.— The X-ray luminosity function (filled circles) for the ROSAT Brightest Cluster 
Sample in the 0.1 — 2.4 keV band. The data are grouped such that each bin contains 17 
clusters, except the highest-luminosity one which contains 12 clusters. The dashed line shows 
the best Schechter function £t. The shaded region marks the luminosity range not covered 
by the BCS. Overlaid is the XLF of the B50 sample of Edge et ah (1990) when converted 
to the 0.1 — 2.4 keV band and grouped such that each bin contains 12 clusters, except the 
highest-luminosity one which contains 10 clusters (open diamonds). The luminosity range 
covered by the B50 in the 0.1 — 2.4 keV band is indicated by the solid vertical lines; the best 
Schechter function fit to the unbinned B50 data is represented by the dotted line. 
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Fig. 2.— The cluster X-ray luminosity function in the 0.5 — 2.0 keV band. The dashed 
curve represents the best Schechter function £t to the BCS data (Filed circles; binning as in 
Figure |l]). The open circles and the dotted line show the group XLF and the corresponding 
best-Ftting power law from BLL. The solid line represents the best Schechter function Ft 
to the XLF obtained by combining the binned BCS data with the Frst 4 data points of the 
groups XLF as determined by BLL. The dot-dashed line, Fnally, shows the XLF determined 
by De Grandi (1996) for a sample of 111 clusters selected from ROSAT data in the southern 
hemisphere. The luminosity range of this sample is indicated by the solid vertical lines (De 
Grandi, private communication). 
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Fig. 3.— The cluster X-ray luminosity function in the 0.3 — 3.5 keV band. The dashed 
curve represents the best Schechter function £t to the BCS data (hlled circles; binning as 
in Figure The diamonds in light and dark shading show the EMSS XLF of Henry et al. 
(1992) in the 0.14 < 2 ; < 0.2 and the 0.2 < 2 ; < 0.3 redshift shell, respectively, while the 
EMSS XLF in the 0.3 < < 0.6 shell is represented by open diamonds. 
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Fig. 4.— The cluster X-ray luminosity function in the 2 — 10 keV band. The dashed 
curve represents the best Schechter function fit to the BCS data (filled circles; binning as in 
Figure]^). The two dotted lines show the power-law representations of the XLF of Piccinotti 
et ah (1982) with and without the Virgo cluster, respectively. The dot-dashed line, hnally, 
represents our best Schechter function fit to the B50 data of Edge et ah (1990) which covers 
the luminosity range bounded by the solid vertical lines. 






Fig. 5.— Likelihood ratio contours (68 and 90% confidence for two interesting parameters, 
i.e., AC = 2.30,4.61) in A and L*x for BCS subsamples obtained by splitting the sample at 
an intermediate redshift ;2sep- The various values of ;2gep are shown in the lower left corner 
of each plot as are the numbers of clusters in the low- and the high-redshift subsample, 
respectively. Filled contours correspond to the low-redshift subsamples; the contours for the 
high-redshift subsamples are shown as bold solid lines. In the Schechter function fits, a was 
kept frozen at its global best-fit value of 1.85. 



